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We consider a quantum field on a lukewarm black hole spacetime. We introduce a new uniform 
approximation to the radial equation, constructed using an extension of Green-Liouville asymp- 
totics. We then use this new approximation to construct the renormalized vacuum polarization in 
the Hartle-Hawking vacuum. Previous calculations of the vacuum polarization rely on the WKB 
approximation to the solutions of the radial equation, however the nonuniformity of the WKB ap- 
proximations obscures the results of these calculations near both horizons. The use of our new 
approximation eliminates these obscurities, enabling us to obtain explicitly finite and easily calcu- 
lable values of the vacuum polarization on the two horizons. 
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I. INTRODUCTION 



The properties of the quantum state in the exterior re- 
gion of a nonextremal static, asymptotically flat, spheri- 
cally symmetric black hole in thermal equilibrium at its 
Hawking temperature is by now well understood ^{5]. 
However, fundamental new issues arise in the investiga- 
tion of space-times with multiple horizons having equal or 
unequal surface gravity /Hawking temperatures or indeed 
multiple coincident horizons with zero Hawking temper- 
ature. Multiple horizons arise, for example, as cosmolog- 
ical horizons in spacetimes with a nonzero cosmological 
constant or as inner horizons in space-times with charge 
or spin. In this paper we consider a scalar quantum field 
on the exterior region of a Reissner-N6rdstrom-de Sitter 
black hole spacetime, where the event and cosmological 
horizon have the same temperature, a so-called lukewarm 
black hole [T]. 

One of the simplest measures of the quantum state of 
a scalar field is the expectation value of the square of 
the field, (<^^), the vacuum polarization, which is of di- 
rect physical significance to the theory of spontaneous 
symmetry breaking ( see e.g., [5]). A central difficulty in 
analyzing {ip'^) arises from ultraviolet divergences which 
need to be renormalized to obtain a physically meaningful 
quantity. In this paper we present a new approach to the 
calculation of {(fP')ren for a quantum field on a lukewarm 
black-hole spacetime in its Hartle-Hawking state repre- 
senting the field in thermal equilibrium at its Hawking 
temperature. (The calculation of renormalized expecta- 
tion value of the stress-energy tensor, (Tpi/)r-en, will be 
presented elsewhere f51.) We note that once the renor- 
malized expectation values in a particular state have been 
calculated, the corresponding results in another state are 
easily computed, as the difference between expectation 



values in two states does not require renormalization 
,as the ultraviolet divergences are geometrical and state- 
independent. 

Previous calculations of {ip'^)ren for lukewarm black 
holes [TIT make use of the WKB approximation. This 
method encounters significant difficulties near the hori- 
zons of the black hole, as will be outlined in Sec. [H] We 
introduce a new uniform approximation, constructed us- 
ing an extension of Green-Liouville asymptotics due to 
Olver [TT], which over comes these difficulties and en- 
ables us to construct {(p^)ren for lukewarm black holes in 
a more satisfactory manner. 

Other approximations have been developed which are 
uniform in I near the horizons [12] . While these approx- 
imations can be used to prove the regularity of {(fP')ren 
on the horizons of a lukewarm black hole [lOj , they suffer 
from two major drawbacks. First, they fail to give exact 
horizon values for {(fP')ren for all parameter sets. Sec- 
ond, and perhaps most significantly, they cannot be used 
to investigate the regularity of (T^i/)ren on the black- 
hole horizons. Our new approximation overcomes these 
shortcomings which, in this paper, will allow us to obtain 
explicitly finite and easily calculable values of {((P')ren on 
the two horizons and in a future paper [^1 will provide a 
route for us to extend this analysis to (T^jy)ren- 

Frolov et al. |13j have developed an approximation for 
2-D static black-hole spacetimes which is valid in the en- 
tire exterior region. However, our approximation is con- 
structed in a systematic manner which affords us greater 
control over the error of the approximation, a property 
which is key for the extension of our analysis to (7)ii/)ren- 

This paper is organized as follows. In Sec. |ll] we will 
outline the existing results and the issues they encounter. 
In Sec, mil we introduce Olver's extended Green-Liouville 
(EGL) asymptotics and demonstrate its application to 
this particular calculation. In Sees. |IV[ |Vj and |VI| we 
describe the calculations needed for {(p^)rem including a 
new method for obtaining explicitly finite horizon values. 
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Sec. VII contains plots of our results. Finally our conclu- 
sions are presented in Sec. |VIII[ Throughout this paper 
we use units in which SttG ~ h = c = ks = ^- 
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II. 
A. 



PREVIOUS RESULTS 

Construction of {(p^)ren 



We will briefly outline the standard method used for 
the calculation of {f(P')ren\ a more detailed description 
can be found in [3]. We begin by considering a scalar 
field (p with mass m and coupling ^ to the Ricci scalar 
R. The field satisfies 



(□ - - iR)(p = 0. 



(2.1) 



We will perform our calculations on a static, spheri- 
cally symmetric background spacetime with line element: 



ds^ = -j(r)At^ 



1 



■^vt?QA(\?. (2.2) 



The spacetime will have a horizon at r = tq whenever 
/(ro) = and in this case the surface gravity of that 
horizon is given by — |/'(ro)|/2 . 

Following the standard procedure we Euclideanize our 
spacetime, that is we perform a Wick rotation r — > ii, 
Eq. (2.2) then becomes 



ds^ = /(r)dT2 



1 



■dr^ 



Assuming kq Oj this space will have a conical singu- 
larity whenever /(rp) = which may be removed by 
making r periodic with period in/nQ. This periodicity 
in the Euclidean section corresponds in quantum field 
theory on the Lorentzian section to a thermal state with 
temperature T — k,q/2-k. 

The unrenormalized vacuum polarization can now be 
written as the coincidence limit of the two point Eu- 
clidean Green's function, which satisfies 



iU^-m' -^R)Ge{x,x') 
so we have 



y/gix) ' 



lim Ge{x,x'). 



(2.4) 



(2.5) 



We are taking our field to be in a Hartle-Hawking state, 
corresponding to the Eucidean solution with the conical 
singularity at tq removed, as such Ge{xtx') takes the 
form j3| 



Ge{x,x')^^ J2 e 



inKQir-r'} 



n— — oo 



xy(2l + l)Pi{,cos-i)CniPni{r<)qni{r>), (2.6) 



/=0 



where Pi is the Legendre polynomial and cos 7 = 
cos cos 0' -f sin6'sin0'cos((/3 — ip'). The functions pni 



and qni satisfy the differential equation 





df\ dS 




dr J dr 






I 





£.R]S = 0, 



(2.7) 



with appropriate boundary conditions, which will be dis- 
cussed in Sec. |Vj This equation cannot, in general, be 
solved analytically and so Pni and q^ji have to be ob- 
tained numerically. The normalization constant C„; is 
fixed by the Wronskian condition 



Gnl 



Pnl- 



dqnl 



qni- 



dpnl 



1 

727- 



(2.8) 



dr dr 

We could, of course, absorb this into the definition of the 
functions but it is frequently convenient to normalise pni 
and qni by their behaviour near different singular points 
of Eq. (pjl). 



The standard procedure to renormalize is to first regu- 
larise by separating x and a;', choosing a temporal split- 
ting as it simplifies numerical calculations. Correspond- 
ingly we set r = r' , 9 = 6' , and (f = cp' and we define 
e — T — t' . Inserting these conditions into Eq. (2.6) we 
have 



r-^de^ + r^siT?eA4)'^. (2.3) Ge{x,x') 



T 
in 



E 



00 

1=0 



{2l + l)GniPnl{r)q„i{r). 

(2.9) 

Due to the distributional nature of the sum, this ap- 
proach carries with it a superficial divergence in the 
Greens's function (2.9), manifested through the noncon- 



vergence of the sum over I in Eq. (2.9). This divergence 
is nonphysical and occurs even when the points are sep- 
arated. Fortunately, due to the high degree of symmetry 
here, it is possible to resolve this issue by subtracting 
multiples of the delta function (which vanish when the 
points are separated) inside the sum, which render the 
sum convergent. In more general space-times, a more 
careful approach avoids these superficial divergences [13] . 
In the present case, the exact form of the subtraction 
term is well known and is given by [3] 

GEix,x') = 

00 00 

^ cos(nKoe) ^ 



1^0 
47r2 



n=l 



1=0 



1=0 



{21 + l)CniPniir<)qujiir>) 
1 



{21 + l)GoiPoi{r<)qoiir>) - 



7^ 



1 

rf\ 
(2.10) 



where for later convenience we have separated out the 
n — term. To renormalize this expression we subtract 
the Christenscn renormalization counterterms and then 
take the limit e — > 0. These terms are given by 

' ' ' V(e-i)i?)ln 



SttV 167r2 ^ 
1 



e2>2|^| 



167r2 967r2^"^ a 



(2.11) 
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Here a is equal to one half the square of the distance be- 
tween the separated points along the geodesic connecting 
them, cr" = CT'", 7 is Euler's constant and is the 
Ricci tensor. For a massive scalar field the constant /i 
is conventionally set equal to m; however, for a massless 
scalar field the constant ji is arbitrary |3]. 
For temporal splitting we have 



1 



1 



while = a'^ = 0. {(p^)div then simplifies to 



/ 



12 



f" 



f 



167r2 1927r2/ 967r2 ASn^r' 



e^V/e2 
4 

(2.12) 



to CniPni{r)quii{r). The idea is that one subtracts this 
approximation inside the sum over which causes the 
sums over I and n to converge rapidly. Then one performs 
the sum of the approximation explicitly and adds this 
back onto the rapidly convergent sum. Of course, this 
does not affect the final answer as all one is doing is 
subtracting the approximation and then adding it back 
on again. 

The standard approximation that is used is the WKB 
approximation and, indeed, previous calculations for the 
lukewarm case adopt this approach [TU]. The WKB ap- 
proximation, however, suffers from problems with uni- 
formity: the closer to the horizon one gets the higher 
I one requires. Thus, if one considers regions bounded 
away from the horizons of the space-time under consider- 
ation then the WKB approximation is useful; however, it 
fails to capture the correct behavior as these horizons are 
approached. This nonuniformity is of particular impor- 
tance when one considers the case that both {ifi'^) numeric 

and {(p )analytic 



We may now express Eq. (2.121 in terms of mode sums us- tion of Eqs. (2.131 and (2.14) and the definition of the 



ing the following distributional identities, valid for small 
e and any ko > ^15] , 



diverge on the horizons. From inspec- 



1 



ncos(nKoe) — + 0{e 



4ln(Va = y^^^i!^+0(e^). 



After subtracting Eq. (2.12) from Eq. (2.10) we may take 



the limit e yielding the result [3j 



surface gravity one sees that this is the case whenever 
m? -|- — is nonvanishing . It can be shown an- 
alytically that these divergences do in fact cancel (fO] 
however the numerical cancellation of these two diver- 
gent quantities obscures calculations near the horizons 
and reduces the accuracy of the results obtained. An- 
other consequence is that it is impossible to find explicit 
values for {((P')ren on the horizons. In addition, it is far 
from obvious how to extend the analytic results to deriva- 
tives of the field as required in calculations of (7)ii/)ren- 



numeric 



nalytic 



with 



numeric 



47r2 ^ [ / 2nKo 



E 



(=0 



{21 + l)CniPni{r)qu.i{r) 



r/5 



Kg 

8^ 



E 

1=0 



{2l + l)CoiPoiir)qoiir) 



rf2 



(2.13) 



In Sec. Ill we will show how replacing the WKB ap- 
proximation by the uniform EGL approximationelimi- 
nates these obscurities enabling us to obtain explicitly 
finite and easily calculable values of the vacuum polar- 
ization on the two horizons. 



III. EXTENDED GREEN-LIOUVILLE 
ASYMPTOTICS 

The fundamental problem that faces us is to obtain 
suitable approximate solutions to a second order differ- 
ential equation of the form 



analytic 



l^K4-(C-i)i?)luf^ 



d w 
da; 2 



= {k^F{x) + G{x)) w, 



(3.1) 



/" , /' 
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B. Previous calculations of {(p'^)re 



(2.14) 



The sums over I and n contained in Eq. (2.13) converge 



so slowly as to make their computation impractical. To 
solve this problem one makes use of an approximation 



with a large parameter k. The form of the approximate 
solution of this equation depends on the zeros and sin- 
gularities of F and G in the region under consideration. 
In the case where there are no zeros or singularities, Li- 
ouville 5I63 and Green independently developed a 
uniform approximation to solutions of such an equation. 
This approximation is widely known in the physics liter- 
ature as the WKB approximation in recognition of the 
development of the theory by Wentzel [18^, Kramers |19j . 
and Brillouin [20) . This approximation remains valid un- 
der some loosening of these conditions on F and G; how- 
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ever, significant reanalysis is needed in the case where 
the region under consideration contains a point xq where 
F{x) has a simple pole and G{x) has a double pole. This 
is precisely the case which arises when the radial equa- 
tion [2]7] is transformed into the current form. The anal- 
ysis to find uniform approximations to differential equa- 
tions which possess such transition points was developed 
by Olver and summarised in his book JJj. We will now 
briefly outline this analysis before demonstrating its ap- 
plication for the case under consideration; a more de- 
tailed description can be found in Chapter 12 of Ref [TT] . 



in which ip{x) = — \{n'^ — and is analytic 

at ^ = 0. This equation has, for each value of k and 
each nonnegative integer j, solutions W2j+i,i(fc,0 and 
W^2j+i,2(fc, C) which are repeatedly differentiable in the 
interval (0,/3), and are given by 

s=0 

+ lln+i{ke")Y.^ + (3.6a) 



A. Background Theory 

We begin by considering solutions of equations of the 
form 



d w 
da; 2 



^ {k^F{x) + G{x))w, 



(3.2) 



in which fc is a large parameter, and at a; = xo, say, F{x) 
has a simple pole while (x — a;o)^G(a;) is analytic. We 
first transform our independent and dependent variables 
in the following manner: 



1 /de 



AF{x), 



w 



da; 



-1/2 



w. 



We assume, without loss of generality, that F{x) has the 
same sign as x — xq- Integration then yields 



1/2 _ 



F^^^{x')dx' ix>xo) 



(-0^/2-/ {-F{x')f^dx' {x<xo) (3.3) 



These equations determine a continuous one-to-one cor- 
respondence between the variables x and ^. The trans- 
formed differential equation is given by 



~de 



4$ 



+ i^ix) w, 



(3.4) 



where 



F{x) 



G{x) 1 d^f(a;)i/4 

^ F(a;)i/4 d^^ 

2 



da; 



A^F{x). 



If G{x) has a simple or double pole at a; = xq, then ip{x) 
has the same kind of singularity at ^ = 0. We denote 
the value of ^'^tp{x) at ^ = by T(n^ — 1), and rearrange 



(3.4) in the form 
d'^W 



d^2 



— 1 ipix) 



4^ 452 



(3.5) 



MO 

^2s 



i-1 



^-K^+^{ke/^) J2 ^ + e2,+i,2(fc, C), (3.6b) 



s=0 



where /„ and Kn are the modified nth order Bessel func- 
tions of the first and second kind respectively. The coeffi- 
cient Ao{£^) is a constant which, without loss of generality, 
we take to be 1 while the higher order coefficients As{(,), 
Bs{^) are determined by the recursion relations: 

= nB,i - ^B'M) + J mSsiOdO (3.7) 

The error terms £2^+1,1 (^,0 and €2^,-1-1,2(^,0 associated 
with each solution satisfy the bounds 



(3.8a) 



|e2,+i,i(fc,OI<A„e'/'/«(fce'/') 

^ \T ) — fc2j+i — ' 

\e2j+iAk,0\<Ke^^Kn{ke^') 

P I T J pj^i ' ^ ^ ^ 



de2j+i,i{k,S,) _ n + 1 
2C 



e2j+i,i(fc,0 



X exp 



K,,{ke'^ 

A 



k 



Vo,c(4 Bo) ] — , 



(3.8c) 



de2j+iAk,C) n + l 

< Ke''K^+,{ke")e^v {^^v,.,(,,){e''B,) 



(3.8d) 
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where Va,bif) is the variational operator defined by 



tions for pni and qni respectively: 



VaMif) = / \f'ix)\dx 



and 



1.07 n = 
1 n>0 



B. Application to black hole spacetimes 



It can be shown that Eq. (2.7) can be written in the 



form of Eq. (3.2). Hence we can apply the above theory 



to find a uniform approximation to either pni or qni of 



order j to Eq. (2.7) by truncating the series Eq. (3.6a) or 



Eq. (3.6b) at an appropriate j. In the spirit of Sec. I 



lA 



we perform the following transformation of variables 



1 



dr' 



(3.9) 



where rg denotes the location of the horizon. Eq. (2.7) 
then becomes 



^ ■' ' s=0 



^2s 



+ hj+iAk,£,{r)), 



(3.11a) 



^ ' s=0 



+ ?2i + l,2(fc,C(0)> 



(3.11b) 



where e = e/(^r^/)^/'*. It can be shown [TT] that for fixed 
r and large k 



<^2j + l.: 



0{k 



1,2 



uniformly with respect to r G [''oj^'b]? provided / does 
not vanish in (rg , r?,] . So we can see that the zeroth order 
extended Green-Liouville (EGL) approximation (j = 0) 
is an approximation to 0{l~^) for large I, the first 0{l^^) 
etc. Therefore for large I 

CniPni{r)qni{r) - Poir)Qo{r) - ©(r^), (3.12a) 
CniPni{r)qni{r) - Pi{r)Qi{r) - 0(r5). (3.12b) 



Inspection of Eqs. ( |2.13 1 and (3.12b) shows that subtract- 
ing the EGL approximation to first order will render the 



with 



and 



(fW_ 



4C 



4^2 



l{l + 1) + (m^ + Riro)) r'o 



+ 

1(1 + 1] 



^ + 2.oro)+^(l-n^) 
-N. 



4,/ 
16 



r (m 

- + 
64/ 



-mr)) 



4 
3 

16C 



16 



4^ 



3rf 
16 

N 

-. (3.10) 



It can be easily shown, using the expansions in Appendix 
|A) that V'(0 is O(0 as ^ ^ (= r ^ tq). From Sec. [niA 
we then have the following uniform jth order approxima 




Event 
Horizon 



Cosmological 
Horizon 



FIG. 1: A comparison plot between the exact solution 
Co2Po2?02 and its zeroth order EGL approximation PoQo for 
a massive conformally coupled scalar field with mass m on 
a lukewarm blackhole with M = Q = 0.1m and A = 3m~^. 
The approximation can be seen to capture the behaviour the 
event horizon while it the breaks down near the cosmological 
horizon. 
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sum over I 0{l/l'^) making the sums converge rapidly. 
A comparison plot between the exact solution and its 
zeroth order EGL approximation PqQq can be found in 

Fig.g 

In our derivation of the EGL approximation in 
Sec. |III A| we implicitly assumed that xq was the only 
singular point of our differential equation in the region 
under consideration. However for a black hole space-time 
with multiple horizons, such as our lukewarm black-hole 
configuration, the radial equation with have a second sin- 
gular point, at ri, say. (For the lukewarm black hole, we 
will take rg to be the event horizon of the black hole and 
ri to be the cosmological horizon.) In this case, our orig- 
inal approximation fails to be uniform as we approach 
ri and we need a second EGL approximation based on 
an expansion around ri. The derivation proceeds along 
the same lines as before yielding the following approxi- 
mations for Pni and Qni , respectively: 



s=0 



\3/4 



-X„+i(fc|(r)V2)^ 



hj+iAk,i{r)); (3.13a) 



f,2s 



e(r 



|3/4 



k,+ia{k.i{r)), (3.13b) 



where 



^l{l + l) + {m'' + {i-l)R{ri))rl 
+ .^(^+2.,.,)+i(l-n^) 

= ;(/ + i) + iv 

fi 1 
C = ( / ^=dr' 



and whe re As and satisfy the recursion relations 
Eq. (3.7) with N replaced by N in the definition of -0. 



The same large I behavior applies to ei_2 as in the previ- 
ous case. 

We now have approximations P and Q to the func- 
tions Pni and qu,i respectively that are uniform in I over 
a closed region including the event horizon but not in- 
clusive of the cosmological horizon and vice versa for P 
and Q. Our approach, which will be discussed in more 
detail in Sec. |IV) then is to subtract and add on both ap- 
proximations by using a switching function to alternate 
between them. 



IV. FORMAL CALCULATION OF {(^2)„n 

We now turn to the details of the calculation of {(fP')ren 
for lukewarm black holes. We will first briefly review 
the key features of these spacetimes before describing the 
details of our calculations. 



A. Lukewarm Black holes 

Lukewarm black holes are a special class Riessner- 
N6rdstrom-de Sitter spacetimes with line element given 
by Eq. with 



fir) 



1 



2M Q2 

r 3 ' 



(4.1) 



where M, Q are the mass and charge of the black hole re- 
spectively and A is the (positive) cosmological constant, 
with Q = M. The Ricci tensor is given by 



P'llU 




(4.2) 



giving a constant Ricci scalar R = 4A. For 4M < y^S/A 
we have three distinct horizons, a black-hole event hori- 
zon at r = r/i, an inner Cauchy horizon at r = r_, and a 
cosmological horizon at r — rc, where 



3/A -1 + a/1 + 4MVA/3 




3/A 1 - Jl -4MVA/3 



(4.3a) 
(4.3b) 
(4.3c) 



The 4th root of / is negative and hence nonphysical. The 
Penrose diagram for this spacetime can be found in |21) . 

The surface gravities of event and cosmological hori- 
zons coincide and are given by 



4Ma 



We shall confine attention to the region r e [rh, rd , which 
has a regular Euclidean section with topology S"^ x S'^ 
|21) . We will take our quantum field to be in a thermal 
state at the natural temperature T — k/2tt where we 
have dropped the subscript for notational compactness. 



B. Computational Strategy 



It is known that the divergence on both the event 
and cosmological horizons for {(p^) numeric is entirely con- 
tained in the n = mode [TUj , so the WKB approxima- 
tion only encounters the problems outlined in Sec II B for 
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this mode. As the WKB mode is easier to work with than 
the EGL approximation, our computational strategy is 
to use a combination of both WKB and EGL approxi- 
mations. We make use of EGL for the n = mode and 
WKB for aU other n modes. For (T)ii/) the corresponding 
strategy requires the use of the EGL approximation for 
the n = 0, 1 and 2 modes. 



C. WKB approximation 

We adopt the WKB approach of Casals et. al [22] which 
we will now briefly outline. We begin by defining /3(r) by 

/3(r) = CniPni{r)q^i{r). 

Then /3(r) satifies the nonlinear differential equation 



dr 
where 



d/3 



1/2 



dr 



1 



4/33/2 



0, (4.4) 



We are looking for a large I and/or large ui approxima- 
tion, so we seek to express fi{r) as an expansion in inverse 
powers of x(r). To keep track of orders it is convenient 
to replace x in Eq. (4.4) by x/e where e is an expansion 
parameter which we will finally set equal to 1 at the end 
of the calculation. We then write 

P{r)^ePa{r)+e^pi{r) + ... 

To balance Eq. ( |4.4| to leading order we require 

Mr) ^ 



Inserting this into Eq. (4.4) for /3 and solving formally 



order by order in e we find a recursion relation for the 
/3„(r). On doing so we obtain the following expansion 

2ri 

A I r In'"" 

(4.5) 



Pn[,l ) - ^2n-|-2m+l 







where a recursion relation for the ^n,m('') can be ob- 
tained from rearrangement of the /3„(r) recursion rela- 
tions 22 . It can be seen from Eq. (4.5 ) that, for large x 
(or equivalently large Z or n if / 7^ 0), /3„ ~ X~^"~^- So 
for fixed r and large I we have the same asymptotic be- 
havior as the EGL approximation. Hence if we consider 
the double sum 



00 00 



^ ^(2; + \)\CM{r)q^i(r) - fi^(r) - /3i(r)], (4.6) 



n=l i=a 



then the summand is 0(^^*) for large I. In addition, we 
see that once the sum over I is completed, the resulting 
summand is 0{n~^\ Therefore in order to co mpute the 
mode sums in {(p^) numeric given by Eq. (2.13) it is suffi- 



cient to subtract either the EGL or WKB approximations 
to the first order. 



D. Incorporating the Approximations 



We now wish to incorporate these approximations into 
expressions for {if>^) numeric and {fp^) analytic- We do this 
in the following manner, which will be justified later on. 
We write 



{'P^)numer^c = A E { T ^ 2^ ^""^ ^ " 
n=l '' 



+ E(2^ + 1) {CnWni{r)qni (r) - /3o(r) - /3i(r)) 

1 



(=0 



E 

1=0 



(2Z-(-l)(/3o(r)-f/3i(r)) 



.Jl/2 



+ ^E(2^ + l)(C'ozPoi(r-)<Zoi(r) 
-5(r)(7o(r) -|-7i(r)) - S{r) {%{r) -f 7i(r)) 

OO 

+ ^ E(2^ + 1) ('5W7i(r) + 5(r)7i(r) 



1^0 



analytic 



(4.7) 



Stt' 



1=0 



1 



1=0 



i^l + l)7o(r) 



(2/ + l)7o(r) 
1 



1 

rf-2 



87r2 



(e-i)i?)ln 



/ 



/2 



Ho 
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here S{r) is a switching function, S{r) = 1 — S{r) and 7 
denotes the EGL approximation valid on the cosmolog- 
ical horizon. The form of S{r) used is irrelevant. The 
only condition on it is that it has a value of 1 on the 
event horizon and of on the cosmological horizon. In 
this case we found the form 



S{r) — tanh 



r-Vh 



to be the most convenient. 

The analytical sums of WKB approximation over I 
have already been computed in jlO] so we will omit the 
details and just present the results 



E 



1=0 



(2/ + l)(/3o(r) + /3i(r))- 



7 



2uj 



{m' + - i) i?) + J^iuj, r) + Muj, r), (4.9) 



where jQ{ijj,r) and Ji(a;, r) are integrals which can be 
found in |10) and are easily evaluated numerically. 
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We now turn our attention the the sum 
1 



E 

1=0 

oo 

E 



(2Z + l)7o(r)- 



r/i 



1=0 

where 



i2l + l)A{r)IoikOKoikC) 



,Jl/2 



(4.10) 



(4.11) 



and we have, for convenience, defined C = f^/^ We com- 
pute this sum using the Watson-Sommerfield formula, 
vahd for any function analytic in the right-hand half 
plane: 



F A 



dA 



7^ 



1 
2 

« / iA- - IdA 

1 + 62^^^ V 2 



Applying this to Eq. (4.101 we can express the sum as a 
sum of two integrals 



Mr)- 
with 



2\A{r)h{kxOKo{kxO 



1 



dA 



4A 



1 + e 



2irA 



A(r)h{hxOKo{k,xO'^\ 



kl = X' + N+^. 



12 ' 



fc?x =-A2+iV+ j^. (4.13) 



Evaluation of Xi (r) 



For large z, we have the expansion 



/„(z)i^„(z) = 



2z 



1 - 



- 1 
8z2 



0(z- 



Using this together with Eqs. (4.11) an d (4^13 ) we see 
that the upper limit contribution in Eq. (4.151 is 0(fc;~^) 
and so we arrive at the result 

Ii(r) = - A{r)kl[Io{koOKo{koO + h{koC)Ki{koC)]. 

(4.16) 

with fcg = + 1/12. 



Evaluation of X-zir) 

We now examine the second integral. First we note 
that kix — (/cq^ — A^)^/^ gQgg ^ero at A — |fco|, hence 
we must introduce cuts in the complex A-plane which we 
take to run from A = — fcp to —oo and from A = fco to 
oo. So we consider the contour of integration shown in 
Fig. [2] We define X2a, ^2b and as the contributions 



e 

I 1 



Ca 



Inspecting Eq. (|4.13|) we can see that 2fcdfc = 2AdA 
id so we can 

Jko 



and so we can change Ii to an integral over k as follows 

dA 1 



2kA{r)Io{kC)Ko{kC) - 



dk r/i/2 



dfc. 
(4.14) 



To compute this integral we first introduce a large k cut- 
off ki and then we take the limit fc; — >■ cx) after the inte- 
gration is performed. We have [23] 



2kIo{kv)Koikv)dk = k^[Ioikv)Ko{kv)+h{kv)Ki{kv)], 



and so Eq. (4.14) evaluates to 

A{r)k^[Io{kC)Ko{kC) 



Ii(r) — lim 



h{kC)K,{kC)] 



-<k, 



(4.15) 



FIG. 2: Contour of integration for X2. 

from Ca, Cb and Cc respectively. 

Considering first, we allow our variable to run from 
to /cq ~ also as the expression is real we can drop the 
TZ symbol. Hence the integral becomes simply 



22a — 



4ArfA 



1 + e- 



^ iAir)Io{k,xOKo{hxC)) 



I-ia contains a divergence as we approach the event hori- 
zon, i.e as C — ^ 0. We choose to take this out of the 
integral explicitly as follows 

rka-e A\fl\ 

ha = I ^^{A{r)Ioik,xOKo{hxC) + A^lnC) 







l + e^ 



Ah InC 



feo-£ 



4AdA 



l + e 



27rA 



(4.17) 
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with Ah = A.{rh) = l/{Krf^). Details of the expansions 
used above can be found in Appendix A. 

I2b can easily be shown to be 0(eln(e)) and so does 
not contribute. 

Lastly we consider the integral l2c- We allow our vari- 
able A to run from fco + e to oo. In this region kix is purely 
imaginary, so we define k = ikix and use the following 
identity [23] : 

Io{ix)Ko{ix) = {Jo{x)Yo{x) + iJo{x)^) 

giving 



1 + e 



2tt\ 



Air)JoikC)Yo{kC) 



with k = yJX^ — k^. Proceeding along the lines of the 
l2a calculation yields the following expression for I2 



1 + e2.A (A(r)/o(fc.AC)i^o(fc»AC) + Ah \nC) 
^ T^^i-lA{r)MkOY,{kO + Au\r.C 



Finally we have the result 



of this section so far giving 



' numeric 



47r2 



00 f 00 

eIe 

71=1 ^ 1=0 



{21 + l)CnlPnl{r)qul{r) 



8tt 



-/3o(r)-/3i(r) + Jo{uj,r) + Ji{uj,r) 

00 

i2l + l)CoiPoi{r)qoi{r) 

i=a - 

- S{r){jo{r) + 71 W) - S{r){'jo{r) + 71 (r)) 

00 

+ si + l)['5(07i W + 5(r)7i(r)], (4.19) 



1^0 



(V-' } analyti 



(5(r)r(r)+5(r))f(r)) 



K + (e - i) R) In f 
167r2 ^ ^ ^' ' \ Ako 



72 



+ 7^ + 



167r2 1927r2/ 967r2 A^ir^r ASir^f 



(4.20) 



Both of which are now manifestly finite on the horizons. 
Clearly {(p^) analytic contains expressions that can only 
be evaluated numerically, so the label "analytic" here is 
used for the purpose of comparison of our results with 
those of existing methods [31 [TD] . 

The sums in {'p'^) numeric are now amenable to numeri- 
cal computation, which will be discussed in the next sec- 
tion. 



AXdX 



[21 + l)7o(r) - — ^ 
r/2 



= -^lnC 
12 ^ 



1 + e 
°° 4AdA 



^A{r)Io{hxOK^{k,xO + ArMC 
A(r)7r 



-MkOyo{kO+Ah\nQ 



l-fe2-A^ 2 
- A{r)kl{UkoOKo{ko(: + hikoOKiikoCn, (4.18) 

and similarly for 70 (r) with C and N replaced by C and 
iV, respectively. It can now easily be shown, using the 
expansions in Appendix |X] that near the event horizon 



Stt 



167r" 



(m2 + (e-i)i?) ln(r-r,0+O(l) 



V. NUMERICAL CALCULATIONS 

To perform the mode-sums contained in {'p'^) numeric 
we first need to find the modes themselves. We do this 



by numerical integration of the radial equation (2.7). We 



perform all of our numerical calculations using Mathe- 
matica [23]. Previous calculations [2-^; relied on pro- 
gramming languages such as C, C-) — h, or Fortran. The 
use of a higher level application greatly speeded up the 
formulation of the code required to calculate the quanti- 
ties required to construct {(fP')ren- 

Turning our attention to calculating the modes, the 
radial equation (2.7) has regular singular points at r = r^^ 
and r = Tc, we apply the standard Frobienus method, by 
writing 



which ensures {(p^) analytic finite on the event horizon. 
The same analysis applied to Tlr) shows that {'f'^) analytic 
is also finite on the cosmological horizon. From this 
we can also deduce that the contribution of the sums 

Ez^q(2' + 1)71 (^) and YlbLai"^^ + l)7l('') to {^'^)numeric 

is finite on both the event and cosmological horizons. 
At this stage it is helpful to bring together the results 



j=0 



ai{r - tq) 



(5.1) 



to obtain the indicial equation a ^ zLn/2. S is either p or 
q, and tq denotes the singular point under consideration. 
We chose our modes so that Pni is the solution regular on 
the event horizon, divergent on the cosmological horizon 
and vice versa for qi^i. The standard procedure, when 
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one is dealing with an asymptotically flat space time, is 
to obtain pni by direct integration of Eq. ( 2.7 1 with initial 



conditions given by the value of the series solution 



(5.2) 



i=0 



at a point r — rg near r — and then using the 
Wronskian condition (2.8) to obtain q^^i. However in the 



case under consideration we have a finite outer boundary 
(r = Tc) hence we can also solve directly for q^j by inte- 



grating Eq. (2.7) from a point outside the cosmological 
horizon r = r„, with initial condition 



i/2 



(5.3) 



The bi satisfy a seven-term recursion relation, given in 
Appendix |B] In order for our numerical integration to re- 
tain a high degree of accuracy, we found it advantageous 
to take our initial point away from the vicinity of the 
horizon. We found that a point about (rc — r/i)/10 from 
each horizon to be most convenient. It is known that a 
series solution about a regular singular point of a differ- 
ential equation has a radius of convergence of a least the 
distance to the next regular singular point [5S] . Since the 
nearest regular singular point to r = Tc is the event hori- 
zon, we see that the series (5.3 1 is valid Vr G (r/j, rd ■ This 
is not the case for the series (5.2) however, as the near- 



est regular singular point is given by the Cauchy horizon 
r = Ti and since \rh — ri\ « |rc — r/j|/10 this poses a 
problem. To get around this issue we define a new Jaffe- 
like series solution [551 for pni . We define a new variable 



which pushes the Cauchy horizon out to infinity and en- 
sures that the nearest regular singular point is r = rc. 
We then have a new series solution for pni 



■j+n/2 



(5.4) 



which is now valid Vr e [r^,rc). The coefficients Cj are 
determined by a seven-term recursion relation, given in 
Appendix [B| 

We are now in a position to calculate the mode sums in 



(4.7). The rapid convergence of the modes sums in (4.7) 



rely on the cancelation of very large numbers. Therefore 
it is neccesary to calculate the mode functions with great 
accuracy. To ensure this we set our code to a precision 
of 24 digits. As was predicted in Sees. |IIIB] and [TVB| the 
sums over I and a;, are 0(/~*) and 0{lo^'^), respectively. 
Finally we turn to the sum 

Yl'iLai'^^ Since we know a priori 

that this sum does not contribute to the divergence on 
either horizon and considering the technical difficulties 



which would be involved in computing this sum analyt- 
ically, we instead simply did it numerically. The NSum 
algorithm in Mathematica computes the sum analytically 
up to a certain user-definable cut-off value of I and then 
uses Richardson extrapolation to calculate the numerical 
tail of the sum. With the optimum choice of the cut- 
off value, the resulting sum had an error w .002%, more 
than sufficient for our purposes. 



VI. {<fi'^)ren ON THE BLACK-HOLE HORIZONS. 

In this section, we introduce a new method of explicitly 
calculating a value for {(p'^)ren on both the event and cos- 
mological horizon. We will perform the calculation on the 
event horizon, the method trivially generalizes to the cos- 
mological horizon. We calculate the terms corresponding 
to {tp'^)numeric and {(fi'^) analytic for radial separation in 
turn. 



A. (^^) 

numeric 

Up to this point we have only considered the case when 
the points are separated in time. The final answer for 

separation. 

Pi^,i vanishes for n > and so has a value oq — ^/ y ' 
on the horizon. Hence separating in the radial direction 
and placing one point on the horizon yields the result 



is, however, independent of the choice of point 



From Eq. (5.2) if we let r — rh we see that 

1/. ^ 



' numeric 



oo 



(2/ + 1) 



ffi(r)-Ai/2(r)i^o(fcC) 

(6.1) 



where we recall from Eq. UaU that A{r) = C/V^'^f- 

To proceed any further we need an expression for 
Qoiir)- We will show that the zeroth order EGL approx- 
imation captures the required local behavior of qQi{r) 
needed to calculate {ip'^)ren on the event horizon. It 
does not, indeed cannot, contain any global informa- 
tion. We find an expression for this global contribution, 
which needs to be evaluated numerically for the space- 
time under consideration. In the following we have been 
inspired by the method due to Candelas [5^5 to express 
qoi (r) in terms of a local analytical term plus a numerical 
global term which can be extended to the case of a luke- 
warm black hole, however we present a variation of that 
method, which provides a more natural way to extend to 
higher orders. 

We begin by noting that qoi{r) can be expanded in a 
series solution about the event horizon as 123 



Qoiir) = - 




aiPoiir) 



ln(e)) (6.2) 
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where Vh = l{l + 1) + r\{m? + £,R) and a; depends only 
onl,e = r~ r,„ ^ f{r^) and ^ /"(r,). 

Next we note that the singular zeroth order EGL ap- 
proxhnation, Qo{r) = A{ry^^Ka{k(^) satisfies 



B. (^^) analytic 



For radial point splitting, with a = 2s^, where s is the 
radial geodesic seapration 



d 
dr 



with 



Vi{r) 



dr 



Viir) 



1 

4(2 



/ 
4 



Qo{r) - 0, 



^2 

T67 



.3) 



r(3/' + r/') 



Applying the standard Frobenius theory [57] to 
Eq. ( 6.3 1 we may find a series solution to Qoir) abou t the 
event horizon. This series is in agreement with Eq. (6.2) 



up to 0{e'^\n{e)) up to a possible diflference in the ai 
term. Hence we may write 



'analyUc = ^ ^(2? + l)A'/\r)Ko{kC) 



1=0 



1 



47r2s2 16 
1 



167r2 967r2 



^(m2 + (e-i)i?)ln(^ 



(6.6) 



A great advantage of radial point splitting is that it 
allows one to integrate the line element directly for s. 
Since t ^ t' , 9 — 9' , and — (f' the line element becomes 

ds^ = dr^// and hence 



1 



dr'. 



(6.7) 



qaiir) = Qo{r) + Pwoi{r) + i?;(r), 



(6.4) 



where Ri{r) denotes the remainder terms and near the 
event horizon is 0(e2ln(e)) as e — >■ 0. (Comparing our 
method with that of Candelas [2^ we see that Ri{r) is 
equivalent to the q term.) The constants /3/ are deter- 
mined by the requirement that qoi{r) is regular on the 
outer boundary and are calculated in Appendix [Cj 

It can also be shown that for I > {r — rh)^^^^ the 
contribution of the terms Ri (r) + ('') cuts off expo- 
nentially in /. Therefore we have an expression for the 
Green's function, valid when one of the points is on the 
horizon 



For a Ricci-flat space time, one can perform this inte- 
gral exactly and hence obtain an expression for s every- 
where without recourse to approximation. For a general 
non-Ricci-fiat space time, this is not the case, but we can 
expand the integrand in the near horizon limit. Integrat- 
ing Eq. (6.7) near the horizon we find (see Appendix [A|) 



{^')analyUc = ^ J^i^l + 1) A^^^ (r)ifo (^C) - ^g^^, 
1=0 



167r2 
167r2 



487r2rh. 



0(e) 



(6.8) 



Ge{x, x') = ^ V(2; + 1) (Qo(r) + aoA) + 0(eln(e)). 

(6.5) 

The convergence of the sum over I of {21 -(-1)/?; presents 
no problem as it can be seen numerically that the /3; 
constants behave like for large I. 



Inspection of Eq. (6.5) leads us to conclude that while 



the expansion (6.4) is of sufficient order to facilate the 
calculation of (if/ rem a higher-order approximation is 
required to capture the behavior needed for the corre- 
sponding calculation of {Tf^u)ren as it will involve deriva- 
tives of the Green's function. 



Finally inserting the expression (6.4) into Eq. (6.1) we 
find that on the event horizon 



' numeric 



^ oo 



{2l + l)Pi. 



with e — r—rh and f'{rh) = f^- To evaluate this quantity 
we first express the e — divergences in terms of infinite 
sums over /, then we bring them inside the existing sum 
over / and finally we calculate the resulting convergent 
sum. First, we transform variable from r to rj using 



T] = cosh 



dr' 



expanding this near the horizon and then reexpressing 
{(p'^)div in terms of r] gives 



Sn^rliri-l) 

167r2^ 

R fh ^2 

967r2 487r2rft 167r2 



(^-i)i?)ln 



c^^ ^'^r\{ri — 1) 



+ 0(77-1). 



(6.9) 
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We have [2S] 

OO 



1^0 

OO 

E 

1=0 



T] — 1 



{I + l)x'+i 



= In 



x-l 



77 > 1, 



X > 1. 



We can reexpress the latter, using the transformation 
X — >■ ^yr] — 1 + 1, as 



E 



1 



In 



^a + i)(v^-i + i)'+i "A 



We incorporate these divergent sums into Eq. (6.8) to 
give, on the horizon, 



analytic 



= um ( y ^ 



N 

a + i)(v^+i)'+i 



/ Al/2( \ 

(21 + 1) ^ifo (A; cosh- 1(77)) - Qiiv) 

\ ao 



+ Ch 
(6.10) 



with 



967r2 487r2rft 167r2 



We now in turn write rj in terms of e and denote the 
resulting summand as F{l,e). Eq. (6.10) then becomes 



1=0 



It is now appropriate to expand F(l, e) as 

F{l,e) ^ Foil) + Fi{xi)V~e + AF{xi,e), 
where 



(6.11) 



Foil) = \imFil,e) 

e— ^0 



(2; + l)[?A(; + l) -ln(fc)] 



TV 

/Ti' 



(6.12) 



6a2 



xi xf 



AFix, e) = Fil, e) - Foil) - F,ixi)V~e, 
with 



(6.13) 
(6.14) 



and ip is the polygamma function. Using the definition 
of a Riemann integral we can write: 



xi+i - Xl) 



/ Fiix)dx= lim y2Fiixi)i 

Jo xi + i^x, ^ 

where Xk is arbritary point e [xi^i,xi]. From the defini- 
tion of X we see that 

/•OO OO 

/ Fiix)dx = limV Fi(a;;)%/e 

1=0 



since Xf^ is arbitrary. Returning to Eq. (6.11) we now 
have 

lim J:Zo Fil, e) = j:r=o Foil) + r F, ix)dx 

+ \\^j:Zo^Fixue). (6.16) 

Since Foil) + Fiix)^/e are the first two terms in the 
Taylor series in e of i^(Z, e) we see that AF(x, e) must be 
of order e or higher, in other words 

AF = €Gixi,e) 

for some G which has an e dependence of 0(1) or higher. 
Applying the Riemann sum argument here gives: 

OQ /-OO 

AFix, e) = %/i / Gix, e)dx. 

1=0 

Fix, e), FqH), and ^"1(2^) are all clearly integrable in a; 
and so Gix, e) is, by construction, also integrable. Hence 
we can conclude that 



00 00 „oo 

lim^F(;,e) =^Fo(0+ / Fi(x)rfa 
(=0 1=0 



(6.17) 



We have confirmed this result by a direct numerical cal- 
culation of Fil, e) as e ^ 0^ 

The right-hand side of Eq. (6.17) can be calculated 



analytically ( see Appendix |D for details) to give 
1 



{0 ) analytic 



d^^(^'^+^^) + 12 + 6^ 



+ Nil-lni^lrh)) 



iS In 



'r{x,l + is)' 
r (x, I — 7(5) 



R 
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(6.18) 



on the event horizon. Here ^(x, a) is the generalized Rie- 
mann zcta function, F is the Eulcr gamma function and 
S ^ y^N + l/12. 

This method is trivially extended to the cosmological 
horizon, yielding a similar expression to Eq. (6.18) with 
the following replacements ^ rc, 6 ^ S, a ^ a and 
/3 — ?> /3 where 
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Combining Sees. |VI A| and |VIB we are in a position 
to write down an expression for ((^^)ren which is vaUd on 
both the event and cosmofogical horizons of a lukewarm 
black hole, namely 



1 °° 



1 



1 



12 ^ 6a'4 



x=-l 



fa , rn' 



R{ra) 

967r2 ' 487r2ra ' 167r2 ' 



(6.20) 



with N' = N or N when = r^, or rc respectively and 
likewise for all other primed expressions. 



VII. PLOTS OF RESULTS 




Event 
Horizon 



jr - ru) 
0.20 m 
Cosmological 
Horizon 



FIG. 4: Plot of {Lp^)ren and it components for a massive, con- 
formally coupled scalar field with mass m with the parameter 
set M = Q = 0.24m and A — 3Th~^ . Here we are near the 
extremal limit M = Q = 0.25rn. In this case we see that 

{(p^) numeric dominates. Once more numeric, {fp^) analytic 

and hence {(p^)ren are clearly finite on both horizons 



We will now bring together the results from the pre- 
vious sections. For the plots of {( p^)ren we calculate the 
quantities {(p^) numeric from Eq. ( 4.19 ) and {if^) analytic 
from Eq. (4.20) at 100 grid points between the two hori- 



zons and take the value of )ren on the horizons to be 



given by Eq. (6.20). For these calculations, unless other- 



wise stated, we express the dimensionful quantities (Q, 
M, A and r) in units of the inverse of mass of the field 
771, which has dimension of length. In Figs. [3] and |4] we 



plot ((p2 



for two different sets of parameters. 



0.03 



0.02 1 



0.01 




Event 
Horizon 



Cosmological 
Horizon 



In Fig. [5] we plot the value of {(p^)ren on the cosmo- 
logical horizon of a lukewarm black hole for a massive 
conformally coupled scalar field with mass m with the 
parameter set M = Q = O.lm as a function of A. As 
A — >• 0, the surface gravity (and hence the temperature) 
also tends to zero. In this limit, the cosmological hori- 
zon is pushed out to infinity. One would imagine then 
that the value of {(p'^)ren on this horizon would be well 
approximated by the value of {(p'^)ren in a deSitter space 
time at a temperature set by the black hole. The line 
element for deSitter space time can be written as j29] 



1=1 



where 



77 = —ae 



(7.1) 



(7.2) 



So for temporal splitting the Green's function for this 
configuration is given by |29j 



1 



^ SeCTTZ/ 



f^(| + z.,|-j.;2;l + A772/(4w')), (7.3) 



FIG. 3: Plot of {(p^)ren and it components for a massive con- 
formally coupled scalar field with mass m with the parameter 
set M — Q = 0.1m and A = 3m~^. It can be clearly seen 
that {(p^)numeric, {'fi^) analytic and hcuce {(p^) ren are finite on 
both horizons 



with 



,2 _ 9 



- 12(m2i?-i +^) 



(7.4) 



and is a hypergeometric function. It is now easily 
shown that for a quantum field in a de Sitter space-time 
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^ }ren{DeSitter) 
^ ) r en 



-0.0005 



-0.0010 



-0.0015 - 




FIG. 5: Plot of {(p'^)ren for both the cosmological horizon 
of a lukewarm black hole and for a de Sitter space time as 
a function of the cosmological constant. As expected they 
appear to be in agreement as A — >■ 0. 



at temperature T 



(0' 



167r2 



= K/(27r) 
1 

247r2a2 ' 



m 



167r2 
1 - Inifi'^a^)) 



12 



(7.5) 



Finally we examine the validity of the large field mass 
approximation to {'f>^)ren obtained from the order to~2 
term in the DeWitt-Schwinger expansion [30] : 

1 



^^{Rap^sR"^''^ - RapR"^)] , (7.6) 



where Rap-yS is the Riemann tensor and Rap is the Ricci 
tensor. In a Reissner-Nordstrom-de Sitter spacetime this 
becomes 



(0' 



1 

167r2TO2 



He 



i?2 



250 - 96MQr + A8APr^)/r'^ - |A)] . (7.7) 



To compare this with our numerical results, we give all 
our dimensionful quantities in units of L = -^/S/A. This 
allows us to vary the mass of the field in our results. Our 
numerical calculations then show that Eq. (7.7) is a rea- 



sonable approximation for, but only for, very large values 
of the field mass m (m « 2QL^^ when M — Q = O.IL). 
This is consistent with the calculations of Anderson [5] 
for a Reissner-Nordstrom space time. 



VIII. CONCLUSIONS 



In this paper we have considered a scalar quantum field 
in a Hartle-Hawking state in a space-time with two hori- 



zons at equal temperatures. We introduced a new uni- 
form approximation to C„iPniqni and have used this ap- 
proximation to calculate {'f^)ren on a lukewarm Reissner- 
Nordstrom-de Sitter space time. The uniformity of this 
approximation allows us to calculate {^p^)ren on both the 
event and cosmological horizon dealing always with finite 
quantities. We have also shown that for small values of 
the cosmological constant, {(fP')ren on the cosmological 
horizon approaches {(fP')ren for a de Sitter spacetime at 
temperature T = k/2i:. 

While other uniform approximations [12j have been 
used to prove the regularity of {({P')ren on the the black 
hole horizons [10) . they fail to give explicit horizon val- 
ues for {((P')ren- Also they cannot be used to investigate 
the regularity of (T^!y)ren on the black-hole horizons. Our 
new approximation overcomes these shortcomings which, 
as we have just demonstrated, allows us to obtain explic- 
itly finite and easily calculable values of {'f>^)ren on the 
two horizons and provides a route for us to extend our 
consideration to {T^^)ren [9,- 

Finally we note that throughout this project we been 
able to perform all our numerical calculations using 
Mathematica [53]. This greatly expedites the code de- 
velopment process and enables easy modification to other 
spacetimes. 
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Appendix A: Near Horizon Expansions 

In this Appendix we describe the near horizon expan- 
sions needed for various calculations throughout this pa- 
per. We begin with 



/«/^^+f^^ + 0(.^), 



(Al) 



with ro denoting the position of the horizon, e = r — vq 
and /q = f'{ro), etc. From Eq. (Al) we can deduce 



Air) 



dr' 2 
C(r) _ 2 



r/(r)V2 rX 



0(e). 



We also need to calculate expansions for the radial 
geodesic distance s from rp to r. First we have 



-1/2 



1 



1/2- 



^1/2 



4/^ 



3/2 



.1/2 



0(6 



3/2X 
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and integrating yields 

f/l/2 R-f •'/^ 

Jo D/O 

The combinations appearing in the Christensen subtrac- 
tion terms are then 



1 



/o , fo 



47r2s2 i67r2e 96n^ 



0(6) 



(/^o + 2/^) 
2ro 



2ro 

0(6). 



Appendix B: Recursion Relations 

In this Appendix we derive the recursion relations 
needed for the initial conditions described in Sec. Ivl 



Firstly we rewrite Eq. (|2.7|) in the form 
d 



dr \ dr J 

- {n^K^ + r^f{l{l + 1) + {rr? + ^Uy )) S = Q. (Bl) 
We define the following expansions 

r'^f = "^(^c - r){r - rh){r - r„)(r - r^), 



dr 



n^K^ + r^f{l{l + 1) + (m^ + ^^Ry) = ^ Vj{r^ - r)^ 

where Aj, Bj an d V,- are easily determined. Then in- 
serting the series (5.3) into Eq. (Bl| and we obtain the 



following seven-term recursion relation for the bi coeffi- 
cients: 



6, = -(^(z + 2 + n/2 - + 1 + n/2 - j)A,6,_ 



J+2 



3=3 



+ + l + n/2 - j)Bj6,_j+i - J2 V,h^,)/G,il,n), 

J=2 j=l 

with G,{l,n) = {i + n/2){{i + n/2-l)A2 + Bl - Vq) and 
we normalise by taking bo — 1/ ^ nr^. 

We now turn to calculating the recursion relation 
needed to determine the coefficients in our Jaffe-like so- 



lution for -pni- We first change variable in Eq. (2.71 to 



X = r — Ti. We then perform the following transforma- 
tion: 



x — >■ u 



X 



and Uix) ^ xl''^{l - ufY{u), 



where Xh = rh — ri. Eq. (2.7) then takes the form: 



uF(u)^ ( uF{u)-^{l - ufY{u) 
du \ du 

-Viniu)il~u)^Y{u) ^0. (B2) 

where 

F{u) = (a; - Xn){xc - x) ^ (-^ x,i) (xc - ) 

VI — u /V 1 — u/ 

and 

+ uF{u) + 1) + (m^ + ^R) + r,) 

with Xn and Xc denoting the position of the corresponding 
horizons in the variable x. 

We now try a series solution of the form 



/2 



1=0 



We define the following expansions: 



{uF)'{l-uf^Y.^^' 



J=2 



(uF) — {uF)(l - uf - 12(1 - u)^(uF)'^ = Vb,u^ 
du ^ 

, 6 

30(ui^)2(l - uf ~ 6{uF)-^{uF){l - uY" = 



3 = 1 



3=0 



where the coefficients Aj , Bj , Cj , and Vj are easily deter- 
mined. Inserting the above expansions into Eq. (B2), we 



obtain the following seven-term recurrence relation for 
the series coefficients q: 

8 

c, = -j + n/2 + 2){i- j + n/2 + l)AjCi^j+2 

3=3 

7 



^(i-.?+7l/2 + l)4Q-J + l 
J=2 

6 6 

~ ^3^1-3)/ Gt{n, I) 



with Gi{l,n) = {i + n/2){{i + n/2-l)A2 + Bi)-Va, and 
we normalise by taking cq = 1/ -^/ktJ. 
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Appendix C: Calculation of Pi 

To calculate the constants /3; we make use of the Wron- 
skian of Eq. ([2j| 



Pi{r)'^qi{r) - qi{r)^pi{r) = i-. 

ar ar j 



We will calculate /3; for the cosmological horizon, /3;; the 
calculation for the event horizon follows identical lines. 
Dividing the Wronskian by qi{r)'^ and integrating yields 
an integral representation for pi (rj) 



Pi{r) = qi{r) 



dr' 



r'y{r')qi{r'Y 



(CI) 



with r/i being the location of the inner boundary. 

Near the cosmological horizon we have a series solution 
for qi{r) 



qi{r) =a + bi{rc-r) + ©((r^ - r)^) 



(C2) 



where a = l/y^trf and the 6; are constants. Hence we 
can rewrite Eq. (CI I in the form 



Pi{r) qiir) 



1 1 
, r'2/(r') V'ZiW " ^ 
qi{r) /"" dr' 

r'2/(r')' 



(C3) 



Now the first integral converges as r ^ while the sec- 
ond integral can be computed exactly as 



We may now obtain an asymptotic expansion for pi (r) 
as r — )■ using ( C3 1 , 



Piir) 



1 



1 



ln(rc - rh) + 



ln(rc - r) 



1 



ln(rc - r„) 



+ 0((r,-r)ln(re-r)). (C4) 



Alternatively, we may expand the cosmological horizon 
equivalent of relation (6.4) 



PQi{r) = P^ir) + Ago;(^) + -R; 
about r = to obtain [2H1 



^(-ln(fc) + ;^ln(2.jr2)-iln(r 



-7 + /3; +0((re-r)ln(r,-r)). 



(C5) 



Comparison of Eqs. (C4) and (C5| provides us with an 



expression for the /?/ which is amenable to numerical cal- 
culation: 



2 M^c - n) + ^ 



1 



ln(rc - rh) 



In (re - r„) 



+ ln(fc) + 7- -ln(2Kr^). 



dr _ 1 



ln(r - r^) + 



1 



1 



2'tn''ri 



ln(r - r„) 



2K»r2 
1 



ln(r - ri) 



ln(rc - r) = /(r). 



Inspection of this expression reveals that we pick up 
a divergence when either of the endpoints of integration 
coincides with a horizon. This gives us the divergence we 
want on the cosmological horizon. However, we also pick 
up an unwanted divergence at — rh] this divergence 
is not real as it is cancelled by a similar divergence in 
the first integral. To show this cancellation explicitly we 
integrate Eq. (C3| by parts to obtain: 



dr' 



1 



1 



, r'''f{r')\qi{r'Y 
I{r') I{r') 



1 



dr' 



„I2 



fir') 



qi{r')^ 
I{r) 



^+1 ^-^I(r')dr' 



^ '^I(r')dr 



/(r') 



Appendix D: Calculations needed for {'p'^) analytic 

Here we present the details of the calculation of 
-^o(0 .fn°° F ijx) dx with Fq and Fi defined by 
Eqs. 0A2\ and |6l3| . 

Firstly we consider the sum, that is 



OO f 

^ (2/ + l)[-ln(fc)+^(; + l)] 



N 

TTi 



(Dl) 



It is convenient to introduce a large I cut off L, evaluate 
each individual sum separately and show that the terms 
in the individual sums which are divergent as L — >■ c» 
cancel, leaving a finite analytic remainder. We start by 
considering the sum 



J2i2l + l)Hk) 



(D2) 



1=0 



To do this we begin by writing fc^ as 



where we have used I{r)/q'^{r) as r — )■ r/j. 
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where 5 = ^1/12 + iV, then we can rewrite the sum iD2h 
as 



+ 1/2) [in ((/ + i) + iS) + In {{l + i) - i5) 

1=0 



We now consider ^ + 1 + Q') — C(-2^i ?) where ( is the 
generahsed Riemann Zeta Function. We can reexpress 
this in the form [28] 



To evaluate the second term in Eq. (Dl ) we consider 



1=0 



Multiplying out gives 



.21 



1=0 



1=0 



g(L + l + g + fc)^~g(gT 



{q + kY 

OO ^ oo ^ L ^ 



fe=0 



— z ln{q+k) 



1=1 

(D6) 

Now we have the identities [28l 



1/^(0 - V-C^ + 1) = y and V'(l) = -7- 



Differentiating with respect to z gives: 



Inserting these into Eq. ( D6 ) and rearranging gives 



dz 



{({z, L + l + q)- Ciz, q)) =J2{q + k)-^ \n{q + k) 



k=0 



and so setting z = — 1 we have 
d 



L 1=1 



dz 



(C(z,L + l + g)-C(z,g)) 



while setting z — Q gives 

d 
dz 



^^{q + k)\n{q + k), 
2=-l k=0 

(D3) 



Now we expand the right-hand side about z — 1 and 
simplify to give 



z + 1 



2(i/)(i + 2) - 1) 



(C(z,i + l + g)-C(^,g)) 



z=0 



J2 Hq + k). (D4) - {{L + 1) - (2L + 1)^(L + 2)) {z-l)+0 {{z - 1)^) 



fe=0 



Combining Eqs. (D3| and (D4), and using the iden- 
tity |28] 



so differentiating this expression with respect to z and 
then setting z = 1 we find 



dx 

yields the result 



C(x,<5) 



= ln(r(<5))--ln(27r) 



^(2Z -I- 1)^'^ + 1) 



(=0 



'-{L + l){L + 2) + {L+lf^j{L + 2). (D7) 



^(2/ + l)ln(fc) = — (C {x,iL+l+5) ~C{x,\ + i5) 



1=0 



+ C{x,L+l-iS) -C{x,^- iS)^ 



The final sum may be performed immediately 



N 



x=-l 



+ i6 



In 



In ' ' ^ ^ 



1=0 



(D8) 



(D5) 



Combining Eqs. (D5|, (D7|, (D8) gives the result 
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J^F^il) ^ ^^{C i^'L + I + iS) + C {x,L + I ~ id)) 



1=0 



i(L + 1){L + 2) + {{L + 1)2 + N)^P{L + 2) 



id In 



■£(|±|±|)ViL(c(,.+.,)+c(.,i-.*)) 



iSln 



r{^-^s) 



(D9) 



r 



We now consider the limit as L — >■ cx); for large z [3T] 
d 



da; 



1 ,2 



In z — 



1 ,2 1 



2 4 
L 

+ T-2+0{z-'), 



z In z + In z 



(DIO) 



where Ei is the exponential integral function. Each term 
decays exponentially in a; as x — )■ oo 28J. At the lower 
limit each term diverges but by expanding about x = 
we see that the combination yields a finite limit 
as X — > giving 



Inr(z) = (z- i)lnz-z- ^ ln(27 



, (Dll) 



V'(z) 



In z - 



2-^ 



+ 0(z-3). 



(D12) 



Using Eqs. (DlOl, (Dll), (D12) and the definition of 5 



it is straightforward to show that in the limit as i — > oo 
the L-dependent terms in Eq. (D9| tend to iV + 1/12. 



Now examining the integral term, we have 
Fi{x) dx = 



0^ N 
^Ki{ax) + -^Koiax) + NEi 



ax \ 



Fiix) 



6a 



4 -7Vlog(V2) 



(D13) 
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